Properties of SAS Josephson junctions in SO(5) theory. 
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We derive the qualitative behavior of superconductor-antiferromagnet-superconductor (SAS) 
Josephson junctions described by Zhang's SO (5) theory. The main differences between these junc- 
tions and conventional SIS junctions arise from the non-sinusoidal current-phase relation derived by 
Demler et al. for thin SAS junctions. Using a simple approximation to this non-sinusoidal function, 
the current voltage relation, Shapiro steps, thermal fluctuation effects and the diffraction pattern in 
a magnetic field are obtained. 



Recently there has been considerable attention focused 
on the theory of Zhang |5J], which attempts to unify 
the d-wave superconductivity in doped Mott insulators 
with the antiferromagnetism of undoped "parent" com- 
pounds. In this theory, the antiferromagnetic Neel vector 
and the complex superconducting order parameter form 
the components of a 5-dimensional "superspin" vector. 
The SO (5) rotational symmetry of the superspin is ex- 
plicitly broken by anisotropy which favors antiferromag- 
netism for undoped or lightly-doped compounds and su- 
perconductivity for more heavily doped materials. The 
transition from antiferromagnetism (A) to d-wave super- 
conductivity (S) is analogous to the spin-flop transition of 
a uniaxial Heisenberg antiferromagnet in a parallel mag- 
netic field. 

From a phenomenological perspective, the most strik- 
ing manifestations of this theory arise in situations where 
the SO (5) order parameter can be rotated from S to 
A or vice-versa by some kind of external field. One 
way of imposing such a field is by proximity to a ma- 
terial in which the superspin direction is strongly an- 
chored, such as in the superconductor-antiferromagnet- 
superconductor (SAS) heterostructure proposed by Dem- 
ler et al. ||. Here we study the properties of such a 
heterostructure in the presence of an external circuit 
and magnetic field. Our analysis, based upon the SAS 
Josephson current-phase relation obtained within SO (5) 
theory has the potential to provide a series of quali- 
tative tests for the applicability of this theory to doped 
Mott insulators. 

The device proposed by Demler et al. consists of a thin 
antiferromagnetic layer sandwiched between two strongly 
superconducting regions (see inset in Fig. |I]). If the 
anisotropy stabilizing the antiferromagnetism in the A 
region is weak, the proximity effect of the S regions will 
align the order parameter uniformly in the S direction. 
For the A region, there is a critical thickness d c , below 
which it remains purely superconducting. For thicknesses 
d > d Cl the order parameter in the A region tips toward 
the A direction, with maximum tipping angle at the cen- 
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FIG. 1. The dc current-phase relation of an SAS junction 
of thickness S = d/d c , as obtained in Ref. 2. 

ter of the region, and the current-phase relation of junc- 
tions is similar to that of conventional SIS junctions. 
For d < d c , a tipping transition occurs at a critical 
value of the supercurrent or, equivalently, of the super- 
conducting phase difference across the junction. As the 
phase difference is varied through the critical value <p c , 
the slope of the current-phase relation changes discon- 
tinuously, from linear in the superconducting regime for 
<f> < 4> c to roughly sinusoidal for <f> > <j) c where the order 
parameter in the A layer tips toward the A direction. 

The current-phase relation calculated in Ref. ^ and 
shown in Fig. [I], can be approximated by 

f (ftW/osin^ 0< H <(p c 

W) = \ (i) 

[ J sin^ 4> c < \4>\ < 7r 

where 4> c defines the discontinuity in the derivative of 
the current-phase relation and Jo = <^ c /(7r<5sin(/> c ), while 
5 = d/d c and d c /ir — <^4, the antiferromagnetic coher- 
ence length. It should be emphasized that although cf> c 
is a useful parameter for identifying a variety of physi- 
cal effects in the SAS junction, it can be related to more 
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physical quantities. Indeed one can show that 
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A useful way of comparing theoretical predictions of 
junction behavior to existing or future experimental data, 
is to analyze the behavior of the junction within a cir- 
cuit, the classic system being the Resistively and Capac- 
itatively Shunted Josephson (RCSJ) Junction |J. Here 
the junction is placed in parallel with a capacitor C and 
resistance R and a constant current is passed through the 
circuit. Such a model takes into account the capacitive 
effects of the junction and also single particle tunneling 
which produces a normal state resistance. In the fol- 
lowing, we explore the properties of the SAS junction, 
using the current-phase relation of Eq. in an over- 
damped circuit where capacitive effects are assumed to 
be negligible, and we contrast these properties to those of 
conventional SIS junctions which obey a sinusoidal cur- 
rent phase relationship (SCPR), given by Eq. (g) with 
0c = 0. 

The dynamics of the overdamped SAS junction (ignor- 
ing thermal effects) are governed by the equation 
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where R is the resistance produced by single quasiparticle 
tunneling and I c i c is an externally applied dc current. 
From Eq. (^|) it is straightforward to derive the dc I-V 
characteristic of the SAS junction from the Josephson 
equation (K/2e)dcf)/dt — V(t), where V(t) is the induced 
time dependent voltage: 
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A dc voltage is not induced in the system unless the 
driving current is greater than Ith-, the maximum of the 
Josephson current Ij [0]. However, unlike a junction with 
a SCPR, the SAS system exhibits two distinct conditions 
for this to occur, depending on the junction thickness. 
For a thin junction with cj> c > 7r/2, the minimum current 
required to induce a dc voltage is Ith = -fosin^ c , while 
for a thick junction with (j) c < tt/2, a dc voltage appears 
only for I th = I . 

For an SAS system with </> c > 7r/2, in the region 
Iq sin 4> c < Idc < Iq , the dc voltage is given by 
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FIG. 2. In the inset we show the dc I-V characteristic of 
an SAS junction of thickness S within the RCSJ model and 
assuming a negligible capacitance. The main body shows for 
different barrier thicknesses the difference between I-V curves 
of an SAS junction and a normal junction. 

where I = Idc/Ioi while in an SAS system with I^c > ^o, 
irrespective of <fi c the analogous result is 
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There are several comments to be made on the two 
expressions above. For a junction with a SCPR, the I- 
V relation has the form V£ c = Rio \JP — \ ||, which 
is hnite for Id c > Iq, and shows ohmic behavior as the 
driving current Idc becomes large. For < <j) c << tt/2, 
Eq. reduces to this form, since in this limit the SAS 
current-phase relation is well approximated by the nor- 
mal Josephson expression. On the other hand, there is a 
clear qualitative difference between Eqs. (||) and (JsJ) as 
4> c becomes larger. For <f> c > tt/2, as Id c ^osin0 c , Vdc 
tends to zero logarithmically. This implies that as the 
driving current is increased from zero, there should be a 
very rapid increase in the induced dc voltage as soon as 
Idc rises above Iq sin0 c . 

By plotting Vdc/hhR against I /Ith, and relating 4> c to 
6 through Eq. (0), we can compare I-V curves for SAS 
systems of various barrier widths. This comparison is 
shown in the inset of Fig. [| Note that as 8 — * 1 (4> c — > 
0), the I-V curve becomes identical to that of a regular 
junction with a SCPR. On the other hand, as the barrier 
width becomes narrower, the dc voltage develops a sharp 
(logarithmic) leading edge and clearly deviates from nor- 
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FIG. 3. Thermal I-V characteristics for an SAS junction 
of thickness S = 0.1 for a range of 7 = I t hh/ekT. 

mal junction behavior. These differences will show up 
most clearly in the differential resistivity. In the body 
of Fig. U we plot the difference between I-V curves of 
the SAS junction and the normal junction, defined as 
V dm = V dc /I th R - V£/I? h R. The maximum deviation 
from normal junction behavior occurs when the driving 
current Id c is just above the threshold current I t h- 

Thermal effects lower the threshold current to zero, 
interpolating between the T = behavior and a linear 
I-V relation as was shown by Ambegaokar and Halperin 
||. The analogous results are shown in Fig. ||for a thin 
SAS junction [Q. Comparison to Ref. Q shows that the 
largest differences between SAS and conventional junc- 
tions occur at low temperatures and small voltages. 

The addition of an ac component to the dc driving cur- 
rent produces Shapiro steps || in the I-V characteristics 
of RCSJ systems at certain values of the dc voltage ||. 
These steps of constant dc voltage for a range of values 
of the dc driving current can be understood as arising 
from phase-locking between the fundamental Josephson 
frequency of the junction or one of its harmonics and a 
harmonic of the applied ac current. This leads us to con- 
sider the addition of an ac driving current to the SAS 
system described above, since the unusual form of the 
current-phase relation may lead to features within the 
step structure of the dc I-V characteristic. In particu- 
lar, the time dependent current Ij(<fi[i\) induced in the 
SAS junction from the purely dc driving current, should 
have a rich harmonic structure which can respond to the 
harmonics of the additional ac current. 

In order to explore this point, we add to the right hand 
side of Eq. (||) an applied current of the form I ac sin ut. 
By making a change of variables u = ut and defining 
A x = 2eI x R/hu where the subscript x = {dc,ac, J 7 th}, 
the dynamics of our circuit can be described by 



— + A 

du 



A 



dc 



A n 



(7) 



and A ac . Steps are expected to occur when (dtfi/du) = n 
(Vdc = nhu/2e) |J, where n is a non-negative inte- 
ger. This corresponds to the previously mentioned phase- 
locking between the fundamental Josephson frequency 
and a harmonic of the applied ac current. In addition, 
sub-steps should appear when (d(j)/du) = n/m, due to 
phase-locking between harmonics of the Josephson fre- 
quency and the applied ac current frequency. 

From our numerical analysis of Eq. (Q), we show in 
the inset of Fig. ^ a typical stepped I-V characteristic 
for a narrow (5 = 0.5) SAS junction with parameters 
A th = 0.8 and A ac = 6.0. The body of Fig. | shows 
the change in step-heights as the ac driving current am- 
plitude is varied (dashed lines). The zeroth step- height 
is defined as the value of Ad c at which (d<f)/du) becomes 
non-zero, while the first step-height is simply the length 
of the vertical step at (d4>/du) — 1, and so on (see inset). 
When A ac — , the zeroth step-height is simply given 
by Ath = 'ieRIth/fioj, where Ith is the threshold current 
required to induce a dc voltage as defined earlier. In this 
figure we also compare the zeroth and first step-heights 
of the SAS junction to a normal junction with a SCPR 
and the same value of A t h (solid lines). 

There are several distinct features which distinguish 
the SAS Shapiro step behavior from that of a regular 
Josephson system. As can be seen in the figure, there is 
a clear linear behavior in the zeroth step-height for small 
values of A ac for the SAS junction. For a junction with 
a SCPR, the step-height dependence on A ac has a quasi- 
Bessel function behavior where the Bessel function is of 
the same order as the step-height For the SAS 

junction, it would appear that the higher harmonic con- 




which can be solved numerically to obtain (d(j)/du) = 
2eVd c /hu as a function of Ad c , for fixed values of A t h 



0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.1 



FIG. 4. In the inset we show the I-V characteristic of an 
SAS junction with barrier thickness 8 = 0.5, A t h = 0.8 and 
A ac = 6.0. The main body compares the variation in zero 
step-heights produced by an SAS junction (circles) to a nor- 
mal junction with a SCPR (squares) for different values of 
A ac - Also compared are the first step-heights of SAS (trian- 
gles) and normal (diamonds) junctions. 
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tent of I j (Eq.[|), allows phase-locking to occur for all 
values of A ac , in contrast to a regular junction with a 
SCPR, where there are zeroes in the step-heights for cer- 
tain values of A ac . Furthermore, the maxima of the step- 
heights are suppressed for the SO (5) case, and there ap- 
pears to be a slight shift in the periodicity of the step- 
height behavior. It is worth noting too, that from the 
higher harmonic content of the SAS current phase re- 
lation, there are large step-heights (sub-steps) at non- 
integer values of (d<j)/ du). These are enhanced over what 
would be expected in the Shapiro steps produced by a 
regular junction with a SCPR. 

As a last illustration of the properties of SAS junctions, 
we consider the isolated junction in the presence of a 
magnetic field. A magnetic field will produce Fraunhofer 
diffraction peaks for the maximum current through an 
SAS junction as a function of the magnetic flux threading 
through it. For the SAS junction in a magnetic field, we 
consider the system to be lying flat along the z-axis with 
thickness d, a height h along the j/-axis and a depth L 
parallel to the x-axis. Applying a magnetic field H along 
the positive y direction implies that the gauge-invariant 
superconducting phase difference across the junction is 
given by <f)[x] — 2irHzox/§o + 0o, where (j>o is the phase 
difference at one end of the junction (x = 0) and $o is 
the flux quantum JTlJ] . In general, the current, / through 
the junction will be given by Jj(<p[x])dx where the in- 
tegrand is the current density Ij(<fr[x])/L. For a junction 
with (j) c > 7r/2, one can show that the current is given by 
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where $ = HzqL is the magnetic flux threading the junc- 
tion. The maximal current I rnax is obtained by maxi- 
mizing / with respect to 4>o- This leads to the following 
equation for </>o, 
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On the other hand, for a junction with c/) c < 7r/2, I max 
is given by the regular Fraunhofer expression I ma x = 
Iq\ sin(7r<I>/<l>o)|/(7r$/$o) for <\> c < \4>o\ < tt/2, and then 
is determined by Eqs. (||) and (||) for \<j> \ < <p c . 

We have determined the maximal current through the 
SAS junction as a function of magnetic flux by 
systematically solving Eqs. (||) and (@) for increasing 
values of the magnetic flux ratio $/ < I > o ■ In Fig- f| we show 
the variation of Imax/ hh with the magnetic flux ratio 
for a variety of SAS barrier thicknesses. Fig. |jj shows a 
general suppression of the diffraction peaks and a strong 
linear dependence before the first diffraction minimum, 
reminiscent of the Shapiro step-height behavior. 

In summary, we propose that the existence of an SO (5) 
superspin vector can be inferred from distinctive features 
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FIG. 5. Diffraction pattern of an SAS junction with barrier 
thickness 5. The regular Fraunhofer pattern from a normal 
junction is given by the 5=1 curve. 

of an overdamped SAS junction, particularly the rapid 
increase of the voltage in the low T I-V curve, the linear 
dependence of Shapiro step heights as a function of ac 
driving current, and the linear dependence of the critical 
current at low magnetic fields. 
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